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Abstract 

Physical principles and mathematical structure involved in deriving an analyti- 
cal representation of the internal structure of the Sun is discussed. For a two- 
parameter family of a non-linear matter density distribution, the run of mass, 
pressure, temperature, and luminosity throughout the Sun is presented in terms 
of Gauss ' hypergeometric function. The system of differential equations gov- 
erning hydrostatic equilibrium and energy conservation for the spherical Sun is 
proved to be a laboratory for the application of special functions. 

1 Introduction 

The structure of the Sun is determined by conditions of mass conserva- 
tion, momentum conservation, energy conservation, and specific modes of en- 
ergy transport through the Sun. For considerations of its internal structure, 



rotation and magnetic fields can be neglected so that the Sun is spherically 
symmetric. It may come as a surprise that much of what has been observed and 
theorized about the Sun can be accounted for in terms of very basic physical 
laws: Newton's laws of gravity and of motion, the first two laws of thermody- 
namics, Einstein's law of the equivalence of mass and energy, Boyle's law and 
Charles' law of perfect gases, and Heisenberg's uncertainty principle. The out- 
line of the theory of the structure and evolution of the Sun has clearly belonged 
to the first half of the twentieth century and is connected with names like Lane, 
Emdcn, Schwarzschild, Eddington, Chandrasekhar, Hoyle, and Fowler (a brief 
outline of this history is contained in Mathai and Haubold, 1988). In the sec- 
ond half of this century, observation and theory of the Sun were greatly refined, 
partly by new observationed techniques and partly by computer simulations 
of its structure and evolution. By and large, the picture of the structure and 
evolution of the Sun seems to be well understood. However, there is a serious 
discrepancy between the theory of how the Sun shines and the most direct ex- 
perimental test of this theory. This discrepancy is more specifically called 'the 
solar neutrino problem' and refers to the fact that the Sun is a volume source of 
neutrinos, particles produced by thermonuclear reactions in the deep interior of 
the Sun, and that the copious flux of solar neutrinos predicted by theory does 
not match the flux detected by experiment on Earth over the past 25 years, 
connected with names like Davis and Bahcall (Bahcall, 1989). Solar neutrinos 
are the only particles that have the ability to penetrate from the center of the 



Sun to the surface almost without interaction with solar matter, escaping freely 
into space carrying the most direct information about physical processes in the 
deep solar interior. The solar neutrino problem, which shall not be the subject 
of this paper, stimulated further studies of solar models, both by employing 
computing machines and using analytical techniques of mathematics. This was 
also the justification of an effort to reconsider the derivation of analytic solu- 
tions to the system of differential equations of solar structure based on the very 
basic physical laws referenced above. Particularly, the solar neutrino problem 
has been considered to be a reason to pursue more actively obtaining analytic 
formulae yielding a description of the gravitationally stabilized solar fusion re- 
actor and showing that methods of the integration theory of generalized special 
functions applied to solar physics constitute a laboratory for the application of 
these functions (Mathai and Haubold, 1988). 

The nuclear reactions which cause the Sun to evolve are sufficiently slow that 
the Sun may be assumed to pass through a series of equilibrium configurations. 
The model for the internal structure of the current Sun may be thought of as 
representing the Sun at an instant of time. Separating the time dependence 
of the evolution of the Sun from the equations governing its internal structure 
allow to replace the time dependent partial differential equations by four simul- 
taneous, non-linear, ordinary differential equations of the first order. The four 
equations represent the radial gradients of mass, M(r), pressure, P(r), temper- 
ature T(r), and luminosity, L(r). Since there are four equations but more than 



four unknown physical variables, one needs additional constitutive equations, 
before the system can be solved: an equation of state for solar matter, a nuclear 
energy production rate, and an opacity law. The full system of equations must 
be solved subject to at most four boundary conditions at the surface and the 
centre of the Sun. These boundary conditions ensure that the structure of the 
Sun can be calculated from the four differential equations but they do not ensure 
that there is a single unique solution (Chandrasekhar, 1939; Stein, 1966). The 
procedure of numerically integrating the solar structure equations takes advan- 
tage of large electronic computers and makes it possible to include a variety of 
detailed physical effects and to vary parameters at will (Bahcall, 1989; Noels et 
al., 1993). A second procedure to provide a solar model, but whose contents 
can be understood intuitively not resorting to numerical techniques, starts with 
Buckingham's theorem that a system characterized by n physical variables can 
be described by an ensemble of n-r dimensionsless products of variables, where r 
is the number of variables whose dimensionless representations are linearly inde- 
pendent. This approach provides a qualitative explanation of the fundamental 
stellar structure equations through dimensional analysis. This analysis suggests 
that more detailed physical and mathematical theories are essentially theories of 
factors of proportionality. They eventually yield numerical values for these pro- 
portional factors because more physical assumptions have to be made (Bhaskar 
and Nigam, 1991). The third procedure treats the solar structure equations by 
rigorous mathematics leading to the Lane-Emden equation which is a second- 



order non-linear differential equation describing the structure of a polytrope gas 
sphere. However, explicit analytic solutions of the Lane-Emden equation exist 
only for values n=0,l, and 5 of the polytrope of index n, not covering the spe- 
cific physical model for the internal structure of the Sun (Chandrasckhar, 1939; 
Horcdt, 1990). 

Additional to the three procedures of constructing solar models referred to 
above, there is an approach to find solutions of the solar structure equations 
by making a specific assumption for straightforward analytical integration of 
these equations. It is possible to obtain analytic solar models by separating the 
hydrostatic component from the energy-transport component of the structure 
equations. For that purpose, an analytic density distribution, namely, that the 
matter density in the Sun varies non-lincarly from the center to the thought 
surface, where the density goes to zero, must be assumed (Stein, 1966; Mathai 
and Haubold, 1988). Then it is possible to integrate the equations of mass con- 
servation, hydrostatic equilibrium, and energy conservation through the Sun. 
Together with the equation of state of a perfect gas, the run of density, p(r), 
mass, M(r), pressure, P(r), temperature, T(r), and luminosity, L(r), is deter- 
mined and can be derived in the form of analytic formulae. The physics of the 
problem requires only three independent boundary conditions: M{r) —> and 
L(r) — > at radial distance r = 0; T — > T = and p — > p = at the ra- 
dius r — Rq of the gaseous configuration. The requirement that p and T tend 
simultaneously to specific values, in this case zero, is only one condition since 



the point at which this occurs is arbitrary. This ambiquity can be removed, in 
principle, by assigning the total mass. The boundary is required to be at the 
point where M(r = Rq) — M Q and this provides the fourth condition. Hence, 
the central density, pressure, temperature, and total rate of energy generation 
are determined as a function of the Sun's mass and radius. However, by as- 
suming an analytic matter density distribution, the energy-transport equation 
of the system of structure differential equations can be satisfied at only one typ- 
ical point in the Sun. The procedure thus established to construct an analytic 
model of the solar interior allows to determine the factors of proportionality 
which remain to be an open problem in the dimensional analysis. The proce- 
dure also reveales that the run of all physical variables for the solar model can 
be expressed in terms of Gauss' hypergeometric function. (Luke, 1969; Mathai, 
1993) 



2 Matter Density Distribution 

For the integration of the system of differential equations governing the in- 
ternal structure of the Sun one has to make a choice for an unknown function 
that still leaves room for physical justification of this choice. By intuition one 
expects that the mass is an increasing function while density, pressure and tem- 
perature are decreasing functions throughout the Sun towards its surface. Thus 
we make a working hypothesis that the matter density distribution p(r) varies 



with the distance variable r, 
(2.1) p(r)= P r I 
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where 5 and 7 are kept as free parameters to satisfy at a later point that the 
density distribution determines properly the mass, pressure, and temperature 
distribution in the interior of the Sun. Equation (2.1) takes into account that the 
chosen density distribution reflects the central value of the density p(r = 0) = p c 
and satisfies the boundary conditon p(r = R Q ) = 0, where R Q denotes the solar 
radius. Also, equation (2.1) implies that p oc M/R 3 , where the constant of 
proportionality depends only on the radial mass distribution and the radial 
distance. 



3 Distribution of Mass 

If M (r) represents the total mass contained within the radius r, and p(r) is 
the density at r, then 

(3.1) ^W =4H r»p(r). 

dr 

Using the assumed non- linear density distribution in equation (2.1), integration 
of (3.1) throughout the Sun leads to the distribution of mass 

where 2-F1 (•) denotes Gauss' hypergeometric function, containing the parameters 
5 and 7 of the matter density distribution in equation (2.1) (Luke, 1969; Mathai, 



1993). Equation (3.2) satisfies the boundary condition M(r = 0) = and can 

be used to determine the central value of the matter density distribution p c in 

equation (2.1) in terms of the parameters 5 and 7 of the chosen model of the 

Sun. The condition M(r = R & ) = M Q in equation (3.2) reveals that 

r , ^ n 3M (f + lXf + 2)...(|+ 7 ) 

(6 ' 6) Pc "4^R| 7! ' 

if 7 in equation (2.1) is kept as a positive integer and subsequently using the 

specific relation for Gauss' hypergeometric function of argument one (Luke, 

1969; Mathai, 1993), that is 

2 F 1 (a, b; c; 1) = T(c)T(c -a- b)/T(c - a)T(c - b). 

Equation (3.3) can be used to select the appropriate values of the parameters <5 
and 7 specifying the solar model with matter density distribution in equation 
(2.1). 



4 Distribution of Pressure 

If g = GM(r)/r 2 is the gravitational force per unit mass at r due to the 
attraction of the mass interior to r, then 

dPjr) GM(r)p(r) 
[ ' dr r 2 

is the equation of hydrostatic equilibrium of the spherical self- gravitating Sun 
with dP(r)/dr being the pressure gradient. The internal pressure produced 
by the weight of the overlying layers increases towards the centre while the 



gas pressure must increase correspondingly to achieve the balance of forces for 
equilibrium. This increase is obtained by inward increases of both temperature 
and density. The internal pressure needed to balance is the gravitational force 
per unit mass (GM/R 2 ) times the mass per unit area (M/R 2 ) which gives that 
P oc GM 2 /R 4 for any spherical body in hydrostatic self- gravitation, where 
the constant of proportionality is again determined by the radial distribution 
of mass in the Sun and the particular radial distance at which P is measured. 
The constant of proportionality can be determined by integrating equation (4.1) 
throughout the Sun by using equations (2.1) and (3.2) for the density and mass 
distribution, respectively. We obtain 



where 2-F1O) is Gauss' hypergeometric function and (— -f) m = T(— 7 + m)/r(— 7) 
is Pochhammer's symbol that often appears in series expansions for hypergeo- 
metric functions. The solution of equation (4.1) given in equation (4.2) confirms 
the condition P(r = Rq) = and gives the central value of the pressure accord- 
ing to the chosen solar model characterized by 5 and 7 in equation (2.1): 
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5 Temperature Distribution 

The simplest theory of solar structure is that of a polytrope. These solar 
models obey an equation of state of the form P = Kp( n+1 ^ n throughout the gas 
sphere. Since temperature does not explicitely occur in this relation between p 
and P, equations (3.1) and (4.1) may be solved independently of the tempera- 
ture and luminosity gradients. This equation of state leads to the Lane-Emden 
equation for polytropic gas spheres, which is an ordinary differential equation of 
second order, but can be reduced, by suitable transformations of the variables, 
to an equation of the first order (Chandrasckhar, 1939; Horedt, 1990). In the 
Sun the density is so low that at the temperatures involved the solar material 
behaves almost as a perfect gas, having molecular weight p, obeying the perfect 
gas law 



where k is Boltzmann's constant and Na Avogadro's number. Substituting 
p oc M/R? and P oc GM 2 /R 4 in (5.1) reveals the dependence of temperature 
on mass and radius of the Sun T oc pM/R, where the constant of proportionality 
depends on the mass distribution and the radial distance. We obtain the detailed 



temperature distribution throughout the Sun by using equations (2.1) and (4.2) 
to rewrite the equation of state given in (5.1), that is: 
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where 2-F1 is Gauss' hypergeometric function. Equation (5.2) satisfies the bound- 
ary condition T(r = Rq) = and allows to determine the central value of 
temperature of the Sun as a function of 5 and 7 contained in equation (2.1): 
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At this point of the procedure to construct a model for the internal structure 
of the Sun by assuming the matter density distribution and subsequently in- 
tegrating the system of differential equations, two remarks are in place. The 
contributions of the radiation pressure to the total pressure and the radial de- 
pendence of the mean modecular weigth \x have been neglected in equation 
(5.1). 

The total pressure P at any point in the Sun is the sum of the gas pressure 
and the radiation pressure, P = P g + P r , where P g is given in equation (5.1) and 
P r = \aT 4 , where a is a constant. Writing P g — [3P and hence P r = (1 — (3)P, 
it follows that P = aT 4 /3(l — 0). This ratio of radiation pressure to gas pressure 



increases towards the center of the Sun but even there the gas pressure exceeds 
the radiation pressure by more than three orders of magnitude. This justifies 
that the radiation pressure has been neglected in equation (5.1) (Chandrasekhar, 
1939). 

The outward flow of energy in the interior of the Sun is driven by the temper- 
ature gradient and resisted by the opacity of the material. The nuclear energy 
generated within the Sun has to continually replenish that radiated away from 
the surface. This energy generation by nuclear reactions causes the solar chemi- 
cal composition to change and keeps the Sun evolving. Since the gas in the solar 
interior is completely ionised, the mean molecular weight /i in equation (5.1) is 
given by \i — (2X + |Y + ^Z)^ 1 , where X, Y, Z are relative abundances by mass 
of hydrogen, helium, and heavy elements (X + Y + Z = 1). The dependence of 
X, Y, Z on the radial distance variable r, which is governed by kinetic equations, 
can not be determined in the procedure of constructing an analytic solar model 
by assuming a matter density distribution as given in equation (2.1). Thus, the 
mean molecular weigth has to be treated as constant in the following. This as- 
sumption does not reflect the situation in the real Sun because nuclear reactions 
have changed the originally uniform chemical composition throughout the Sun. 



6 Nuclear Energy Generation Rate 

Nuclear energy production in the Sun depends steeply on the temperature 
of the material and is very concentrated towards the center of the Sun. This 
is one of the major reasons that calculations of the internal structure of solar- 
type-stars made already considerable progress before the physical mechanism 
of the production of energy by nuclear reactions was understood. The rate of 
nuclear energy generation can be written (Mathai and Haubold, 1988), 



where eo is a physical constant depending only on the chemical composition of 
the solar material and the units chosen. Substituting p(r) and T{r) in equation 
(6.1) by equations (2.1) and (5.2), respectively, and taking advantage of the 
fact that < (;^) < 1, the energy generation rate in equation (6.1) can be 
represented in the form of a polynomial 



(6.1) 
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Note that the representation of e(r/i? ) in equation (6.2) is essentially deter- 
mined by the four free parameters 5, 7, n and m in equation (2.1) and (6.1) with 
the only restriction that 7 be a positive integer. The coefficients ao, oq, . . . a2 7 
in equation (6.3) are determined by the following polynomial of degree 27 in 
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7 Luminosity Function 

Let be the function representing the flow of integrated radiation across 
a sphere of radius r. If e is the energy produced per unit time by nuclear 
reactions in each unit mass of solar material, than the balance between energy 
generation in the solar interior and energy loss through its surface is governed 
by the equation of energy conservation 

dL(r) 



(7.1) 



dr 



4irr p(r)e(r), 



where e(r) is given by equation (6.2). Integrating equation (7.1) over the Sun's 
interior leads to the luminosity function in terms of Gauss' hypergeometric func- 



tion: 
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s* = s + ^(3 + 2q) + m + 2n 2 + . . . + (27)n 27 , 

where / is given in equation (6.3) with s substituted by s* . Equation (7.2) 
satisfies the condition L(r = 0) =0 and gives for the total energy output 
L(r = R Q ) = L Q , 
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8 Conclusions 

Assuming a two-parameter family of matter density distributions in equa- 
tion (2.1) made possible the analytic integration of the differential equations 
for conservation of mass, momentum, and energy throughout the Sun, equa- 
tions (3.1), (4.1) and (7.1), respectively. This procedure shows that hydrostatic 
equilibrium, equation of state, and overall energy conservation determine the 
state of the central solar conditions, particularly the gravitationally stabilized 
solar fusion reactor. The mathematical method chosen reveals the factors of 
proportionality which are kept undetermined in dimensional analysis commonly 
pursued to understand astrophysical relationships between global parameters of 
the Sun. A common mathematical element of the derived distributions of mass 



(3.2) , pressure (4.2), temperature (5.2), and luminosity (7.2) throughout the 
Sun is Gauss' hypergeometric function 2-F1 (■) which is numerically easily acces- 
sible through programms for doing mathematics by computer like Mathematica 
(Wolfram, 1993). 

It has been emphasized above that the assumption of an analytic matter 
density distribution means that the equation for the transport of energy by 
radiation through the Sun can be satisfied at only one specific point in the Sun. 
The flow of radiant energy per unit area through the Sun is proportional to the 
ratio of radiation pressure gradient and opacity per unit volume. That is 

d{\aT±)/dr T 3 dT/dr 

(8.1) H oc — oc , 

up up 

where n denotes the opacity per mass unit at temperature T and density p. 
Because the energy flowing out through the Sun is transported by radiation, we 
find for the luminosity L 

T , t,2tt R 2 T 3 dT/dr 

8.2 L oc 4-kR H oc '- — . 

up 

Since for a given solar structure p oc M/R? and T oc M/R, it follows for L that 

(8.3) Loc-M 3 . 

For solar composition, Kramer's power law approximation for the opacity given 
by 

(8.4) k oc K pT~ 7/2 , 

where k is a physical constant depending on the chemical composition of the 
solar material and the units chosen, which leads to a luminosity -mass- radius 



relation, 



(8.5) 
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The differential equation governing the outward flow of energy driven by the 
temperature gradient and resisted by opacity in (8.1) is 



where c denotes the velocity of light. To satisfy the radiative energy transport 
equation (8.6), taking into account the density distribution assumed in equation 
(2.1), the temperature gradient dT/dr in equation (8.6) has to be equal to the 
temperature gradient in equation (5.2). This condition can be satisfied at only 
one specific point in the solar interior, for example at the boundary of the nuclear 
energy producing core region (at r w O.3i? where L « L & ). 
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